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The Reversible Nearest Particle System 
on a Finite Interval* 

Dayue Chen, Juxin Liu and Fuxi Zhang^ 

Abstract: In this paper we study a one-parameter family of attractive reversible 
nearest particle system on a finite interval. As the length of the interval increases, 
the time that the nearest particle system first hits the empty set increases in different 
order, from logarithmic to exponential, according to the intensity of interaction. In 
particular, at the critical case, the first hitting time increases in a polynomial order. 
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1 Introduction 

A nearest particle system on S = {1, 2, • • • , N} is a continuous time Markov chain 
with the state space {A : A C S}. The jump rates are specified as follows: 

q(A, A \ {x}) = 1 if x G A; 

q(A, A U {x}) = (3(l x (A),r x (A)) if x e S \ A; 
q(A,B) = 0 otherwise. 

Here l x (A) and r x (A ) are the distances from x to the nearest points in A to the 
left and right respectively, with convention that l x (A) (or r x (A)) is oo if y > x (or 
y < x, respectively) for all y G A. We assume that 

1. f3(l, r ) = P(r,l); 

2. (3(1, r) is decreasing in l and in r; 
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3. 0(oo, cx)) = 0 ,/3(l, oo) > 0; 

4. Yli 0(1, oo) < cx). 

There are many choices of 0(-,-) satisfying the above assumptions. 

Example 1. (The 1-dim contact process) 0(1,1) = 2A, 0(1 ,r) = (3 {1,1) = A for 
l,r > 1, and 0(1, r) = 0 otherwise. 

Example 2. (The uniform birth rate) 0(1, r) — \/(l + r — 1). 

Example 3. (The reversible case) 

0(1, r) = A ^y? , 0(1, oo) = 0(oo,l) = \fi(l), 

i>(L + r) 


where 

OO 

^(•) >o, XI 

n=l 

Assume further that 


1, and 


ip(n) 


^(n + 1 ) 


\ 1 - 


( 1 ) 


X n2 V , ( n ) < °°- (2) 

n 

For example, V’( n ) = cn~ a for some a > 3 satisfies the above requirement. 

It is helpful to associate a subset A of S' with an element £ of {0, l} 5 and use 
them interchangeably: £(x) = 1 if and only if x G A. Configuration £ will be 
given an occupancy interpretation. We say there is a particle in x if £(x) = 1, and 
we say the site is vacant if £(x) = 0. Then the above transition mechanisms can 
be interpreted as follows: Each particle disappears at rate 1 independently, and a 
particle is born at vacant site x at rate 0(l x (A),r x (A)). 

The transition mechanisms also make sense if we replace {1, 2, • • • , with the 
integer lattice Z. However, because of Assumption 4, the state space (0,1} Z consists 
of four disjoint parts: 

1) all finite subsets of Z; 

2) all subsets of Z with infinite many particles both to the left and to the right of 
the origin; 

3) all infinite subsets of Z with finite many particles to the right of the origin; and 

4) all infinite subsets of Z with finite many particles to the left of the origin. 

The first two cases are extensively studied, and are called finite and infinite nearest 
particle systems respectively. A comprehensive account can be found in Chapter 7 
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of Liggett (1985). The last two cases share many properties of the first two cases, 
and are indispensable in some occasions, e.g., Lemma 4.1. 

For interacting particle systems people are most concerned with the existence of 
phase transition and the critical value. For the infinite nearest particle system with 
the uniform birth rate (Example 2), the critical value is 1, see Mountford (1992). 
For the reversible nearest particle system (Example 3), the critical value is also 1. 
For the contact process (Example 1), the critical value is unknown but is between 
1.5 and 2 , and is denoted as A c throughout this paper. 

Can the critical value of an infinite model be detected by the counterpart on a 
finite interval? This interplay was first explored for the contact process in a series 
paper by Durrett et al. The main results are summerized as follows. Let {Q Y : t > 0} 
be the contact process on {1, 2, • • • ,N} with the parameter A starting from all sites 
occupied, and r ^ be the first time it hits the empty set. 

Theorem 1.1 (i) If X < X c , then there is a constant 71 (A) > 0 so that as N — > 00 , 
T/v/log N —■> 1 / 7 i(A) in probability (Durrett and Liu (1988), Theorem 1). 

(ii) If A > A c , then there is a constant 72 (A) > 0 so that as N 00 , (log tn)/N — > 
72 (A) in probability (Durrett and Schonmann (1988), Theorem 2). 

(hi) If A = A c and a,b e (0, 00 ), then P (aN < tn < bN 4 ) —> 1 as N —» 00 
(Durrett et al (1989), Theorem 1.6). 


We believe that these statements hold for a large class of interacting particle 
systems. In this paper we like to study the asymptotical behavior of the hitting 
time of the reversible nearest particle systems (Example 3) on a finite interval, 
as the length of interval increases. The results read as follows. Let {Cjv ■ N > 1} 
be any sequence of increasing numbers such that linry^oc Cn = 00 . 


Theorem 1.2 Suppose the initial state is {1, 2, • • • , iV}. 

(1) If 


A < min{l, min 


ip(n) 




E/+r=n<K0# r ) 

then Eg n < ClogN for some constant C which is independent of N, and 


lim P(ctn < CwlogiV) = 1 ; 

AT—>00 

(2) If A > 1, then there is a constant 7 > 0 such that liniAr^ooP [a N > e lN ) = 1. 
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Remarks : It is not difficult to establish estimates of the opposite direction, see 
Theorems IQ and 12.41 Together we have shown that a^ increases logarithmically if 
A is small enough and exponentially if A > 1. 

For any non-empty set A = {.x'i, x 2 ■ ■ - , £&}, we assume without loss of generality 
that x\ < X 2 < ■ ■ ■ < Xk and define 


v<p(A) 


4>(x 2 - xi)iIj(x 3 - x 2 ) ■ ■ -^{xk - x k -i) if k > 1; 

1 if k = 1. 


Let S N = {0, K n = S Ae5jv \ {0} v*i>(A), and n{A) = u^(A)/K N . Then it is 

a probability measure on Sjy- 


Theorem 1.3 Suppose that A = 1 and the initial distribution is n. Then 


lim P 


N 


< & N < CnN~ I — 1 


iV—»+oo \ C tv 

We now proceed to prove Theorems 11.21 and 11.31 by three different approaches. 


2 Comparison by Coupling 

We will prove the first part of Theorem II.21 bv establishing a more general conclusion 
(Theorem 12.211 . Let [X t : t > 0} be a birth and death process on {0,1, • • • , N} with 

death rate : a % — i, for i — 1, ■ • • , N\ 

birth rate : b t = (i + l)a, for i — 0, ■ • • , N — 1. 

Let r = inf{t > 0 : X t = 0} be the first time that {X t : t > 0} hits 0. Let E N be 


the conditional expectation on 

X 0 = N. 



Lemma 2.1 Suppose that A" 0 

= N. For large N, 




' (2 log N )/(1 — a) 

if a 

< 1 

E n t < < 

2 N log N 

if a 

= 1 


a N a/(a — l) 2 

if a 

> 1 


Furthermore 

E n t 2 < 2 (E n t) 2 . (3) 
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Proof. Let P l be the conditional probability distribution on the initial state i, E l 
be the expectation with respect to P\ and m l = E 1 t for i = 0, • • • , N. It is shown 
in Wang (1980) that 

N N 

E N T = ^ e n EN T 2 = 22 £i ’ 

2=1 2=1 

where 

1 

e-i = — 

= -(1 + ol + a 2 + • • • a N ~ l ). (5) 

i 

^ _ 2 rrii ^ 26j6j +1 • • • 6j + fcmj_|_fc + i 

1 ' ' ' Q'i+kQ’i+k+l 

k =U 

Notice that m l <rriN for any i < N. It follows that £* < ‘hn^ei. Therefore, 

N N 

E n t 2 = '22 £ i — 2 rriN 22 e * — ^ m NE N r = 2 ( E n t ) 2 . 

2=1 2=1 

If a — 1, by ©, e* = (N — i + 1 )/i, and for large N, 

N N 

E n t = 22 e i < E 22 rl < 21Vlog N. 

i= 1 i =1 

If a < 1, E n t = Eir ei < (1 - a)' 1 Eli < _1 < (2 log JV)/(1 - «); 

If a > 1, then E n t = Eili e * < (a — 1) _1 Eili a w ~* +1 < a 7V+1 /(a — l) 2 . □ 

Consider a nearest particle system {ff J : t > 0} on {1, 2, • • • , IV} starting from 
{1,2, ■ ■ • ,1V} (not necessarily reversible). Let <7 at be the first hitting time of the 
empty set by , and 

M — max {max Y^ @(l,r), Y^ /3(/,oo)}. 

l-\-r=n l 

Theorem 2.2 Suppose the initial state is {1,2, ,1V}. If M < 1, then Ea n < 

(2 logIV)/(I — M); and for any sequence {Cn : N > 1} of increasing numbers such 
that linpv^oc CV = oo, limjv -,.00 P(<Jn < CwloglV) = 1. 

Proof. Let |H| be the cardinality of set A. For any configuration £ such that |£| = i, 
there are at most i + 1 intervals of consecutive vacant sites, separated by occupied 
sites; the rate that a new particle in each interval is born is no more than M. Hence 


JV-l-i 

£ 

k =0 


^ 2 ^ 2+1 * ' * bi-\-k 
^2^2+1 ’ * * ^2 +k ^2+fc+l 


(4) 
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the rate that \£^\ increases by 1 is no more than (i + 1 )M. On the other hand, 
when |£f| = i, the rate that |£ f | decreases by 1 is equal to i, the total number of 
particles. Compare \ft] with a birth and death process X t with parament a = M. 
Since initially X 0 = (£^1, there is a coupling of {X t : t > 0} and : t > 0} such 
that 

^"(^>1^1, Vt>0) =1, ( 6 ) 

where P N >€ is the coupling measure with the initial state (N,f N ). By 0, Cn is 
stochastically dominated by r, i.e., for any t > 0 , 

P(a N >t)< P N (r > t ). (7) 

By the Chebyshev inequality and ©, for any cn > 0 , 

E n t 2 2 

P(a N > c n E n t) <P n (t > c n E n t) < --—— 5 < —. ( 8 ) 

(cnE^t) c n 

For any sequence Cn —> oo as At —> oo, choose cn = Cn{ 1 — M)/2. Then an upper 
estimate of <7 jv may be taken as CnE n t , and the claims in Theorem 12.21 hold by (JHJ) 
and Lemma 2.1. □ 

By the same argument it is not difficult to establish following estimates, though 

a renormalization argument is used in the proof of the second part of Theorem 12.41 
We will skip the proof, since they are not needed in proving Theorems 11.21 and 11.81 

Theorem 2.3 Suppose the initial state is {1, 2, • • • , N}. 

(1) If M = 1, then Ecjn < IV log N; and linrv^oc P (o 'at < CnNIoqN ) = 1; 

(2) If M > l, then Ectn < M N+1 /(M — l) 2 ; and there is a constant 71 > 0 such 
that 

lirn P (<t n < e llN ) = 1. 

Theorem 2.4 Suppose the initial state is {1, 2, • • • , N}. 

(1) For any e > 0, limjv->oo P{&n > (1 — e) log At) = 1 ; 

(2) //max„min{| 3n — l), 00 )} larger than the critical value 

of the contact process on Z, then there is a constant 7 > 0 such that 

lim P (a N > e lN ) = 1. 

JV—>00 v ' 
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3 A Lower Estimate of a\ 

We first extend the notation introduced before Theorem Q For any non-empty set 
A = {xi,x 2 • • • , Xfc}, X\ < x 2 < ■ ■ ■ < Xk, define 

I X k ~ l ^{x 2 - xi)^(x 3 - x 2 ) ■ ■ ■ if(x k - x*_i) if k > 1; 

"4’A a ) = \ 

[1 if k = 1. 

Let 5 tv = {0, - K n ( A) = £ Ae5jv \{ 0} and = ^,x{ a )/ k n( A). 

Then 7r is a probability measure on dTv- 

Lemma 3.1 Kn(X) > CN 2 e"'^ N for X > 1, where y(l) = 0 and'y(X) > 0 if X > 1. 


Proof. 

[iV/3] JV 

a- n (a)= £ E E * lsl ~vo. o) 

?e5jv\{0} *=0 j/=[2JV/3] Z£S N (x,y) 

where 


SN(x,y ) = {£ G Sv : f(x) = f(y) = l,f(s) = 0, V 1 < z < i, or y < z < N}. 


In light of 0 , by the Renewal Theorem, u 1 p(SN(x,y)) > 1/(2 whenever 

y — x is large enough. If A = 1, then ALy > CN 2 when N is large, and we are done. 
If A > 1, we can choose constant 5 > 0 such that 

W((£ e S N (x,y)]\£\ > 5\y-x\}) > 


This together with © implies the desired conclusion. In particular we may choose 
7 (A) = (5/3) log A. □ 

We now use an idea in proving Theorem 7.1.20 of Liggett (1985) to prove 


lim P n 

N—>oo 



Kn \ 

c n n) 


l. 


( 10 ) 


The first half of Theorem 11,31 readily follows from HUD and Lemma o Notice 
that the hitting time of the nearest particle system starting from {1,2, •• • , N} is 
stochastically larger than that starting from the initial distribution it. Therefore the 
second part of Theorem 11.21 also follows, with a little change in 7 . 
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Proof of (1101) . The reversible nearest particle system Y : t > 0} is a Markov 
process taking values in S^ with jump rate 


q{A,B) = 


y4>(l x (A))ip(r x (A)) 

i!>(lx(A)+r x (A)) 

0 


if x G A, B = A \ {a:}; 
if x ^ A, B = A U {a;}; 
otherwise. 


ft is reversible with respect to n in the sense that 7r(A)q(A, B) = it(B) q(B, A) for 

A,BeS N \ {0}. 

Let '■ t > 0} be a Markov process on Sn, which is a modification of : 
t > 0} so that particles can be born from the empty set. More specifically, the 
transition rates of {^ v : t > 0} is defined as follows. 

' q(A,B) if A + ^ 

q(A, B) = < q if A = 0 and \B\ = 1; 

0 otherwise, 

where q > 0 is a constant to be determined later. Let Kn stand for Ku( A), 

v-ib ({0}) = q~ X , and n = v^j (. K N + q _1 ) . 

Then Yr : t > 0} is reversible with respect to n in the sense that Tr(A)q(A, B) = 
7 t(B)q(B, A) for any A, B e Sn- 

Let P be the distribution of {^ Y : t > 0} with initial distribution ff, and E be 
the expectation with respect to P. Notice that {^ v : t > 0} is stationary under P. 
For any t > 0, 

r-2t 

2t7r({0}) = E 1 { gy =0} ds. 

Introduce the stopping time r = inf{t > 0 : — 0}. By the Strong Markovian 

Property, the right side above equals 


-2 1 


EE 


> EE 1 


1 - (j Q 

{£*=0} 


r»2 1 


rr+t 


L {r<q 


1 — d q 

{&=*} 


{ T< t} i 

r t 


Tr \ >EE ^1 

Tr ] = P(t < t)E ( / 1 




{^=0} 


ds 


yAr _ 
so — 


Denote by cr the first time : t > 0} jumps. Then 


E ( f l^ Ar _g^<is 


tTV _ 

So — 


> E =0J = J s %e q9S ds, 












where q 0 = 5^g(0, £) = Nq. Hence 

P (t <t \< 2 t ^(W) = 2 t K l _^_ 

“ f* s%e~^ s ds K n + q- 1 1 - e~ Nqt - Nqte~ Nqt ' 

On the other hand, 


( 11 ) 


P(r <t) > P(r<t, ffti)=P 7 ^)P{r< t\t$ ff 0 ) 


K 


N 


Kn + q 

This together with m yields that 


—P(<7n < t). 


P (cT N <t) < 


2 tN 


K n (1 - e~ Nqt - Nqte~ Nqt ) 


Let q —> 00 , then 


P {a N <t)< 


2tN 


K 


N 


This implies m , by choosing t = K^/(CnN). 


□ 


4 The Critical Case 

In this section we will prove the second half of Theorem II .Ml he., when A = 1, 

lim P n (a N < C n N 2 ) = 1. (12) 

Let {r)t : t > 0} be an infinite reversible nearest particle system on Z with finite 
many particles to the right of the origin (The third case on page 2); and r t the 
rightmost particle in {rj t : t > 0}, i.e. r t := sup{x : r] t (x) = 1}. The properties of 
r t of the critical nearest particle system are studied in Schinazi (1992). For a recent 
survey, see Mountford (2003). 

Lemma 4.1 (Schinazi (1992), Theorem 1) Let {r] t : t > 0} be the critical reversible 
nearest particle system on Z. Suppose the initial configurations have a particle at 
the origin and no particle to the right of the origin, and follows the renewal measure 
Renfip) with density Then, as a —> 00 , r a i t /a converges in distribution to a 

Brownian motion with diffusion constant D > 0 in the Skorohod space. 
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Proof of (1121) . Partition the configuration space Sn according to the position of 
the rightmost particle. Namely, let 


A x = {£ G Sn : £(#) = 1, and £(y) =0 for any y > x} 


be the set of configurations whose rightmost particle is at x. Denote by P the dis¬ 
tribution of : t > 0} with initial distribution 7r, and by Pn.x the conditional 
distribution of the nearest particle system on {1, 2, • • • , N} whose initial configura¬ 
tions are in A x . Then 

N 

P = Y j P(A x )P n , x . (13) 

x=0 

Denote by P the distribution of the nearest particle system on Z with the ini¬ 
tial distribution in Lemma rm and P x the translation of P by x. Thanks to the 
attractive property, there is a coupling of P x and Pn, x such that for all t > 0 and 
all i G Z, 

&(}) < *(<)■ (14) 

Then under this coupling, ^ = 0 once r t < 1, hence < inf{t : r t < 1}. 

Suppose that linpv^oo Cn = oo. For any C > 0 and large N, 


Pn,x (&N < CnN 2 ) 


> Pn,x ( cr N < C(x - l) 2 ) 

> P x (3 t < C(x — l) 2 s.t. r t < l) 

= P (3 t < C(x — l) 2 s.t. r t < —(x — 1)) 
= P (3 t < C s.t. r (x _ i) 2 t /(x - 1) < -l) . 


Here the first equality holds because P x is the translation of P by x. This together 
with Lemma 14.11 implies that 

liminf P Nx ( <j n < C n N 2 ) > P (3 t < C s.t. B t < —1), V C > 0, 

N,x->-\-o o ’ 

where {B t : t > 0} is a Brownian motion with diffusion constant D. Let C —* +cx), 
the right side of the above equation converges to 1. Hence 

lim P Nx ( a N < C n N 2 ) = 1. 

N,x —>+oo x 

Consequently, for any £ > 0, there exists N 0 > 0 such that for any N > x > N 0 

Pn,x (<?N < CjyN 2 ') > 1 — £. 
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This together with m implies that 


N N 

P (a N < C n N 2 ) = {°n < C n N 2 ) > (1 - e) ^ p ( A x)- (15) 

x=N 0 


x=l 


On the other hand, 


No — 1 No —1 x 

^(S N (y,x)) < iV 0 2 . 

x—\ x=l y= 1 


Therefore, as N —>■ oo, 


N 

P( A x) > 1 — Nq/(CN 2 ) —> 1. 

x=N 0 

This together with (USD implies that lini inf lV ^oc P (&N < CnN 2 ) > 1 — e. Let e —> 0 
and the result follows. □ 
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